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We investigate two points related to existing treatments of isospin-breaking cor- 
ps) \ rections to the CVC relation between a[e + e~ — > 7r + 7r~] and oT{t~ — > i/ T 7r _ 7r ]/(is. 

Implications for the value of the hadronic contribution to = {g — 2)^/2 based on 
■ those analyses incorporating hadronic r decay data are also considered. We conclude 

that the uncertainty on the isospin-breaking correction which must be applied to the 
^ t decay data should be significantly increased, and that the central value of the p-uj 

0^ \ "mixing" contribution to this correction may be significantly smaller than indicated 

by the present standard determination. Such a shift would contribute to reducing 
the discrepancy between the r- and electroproduction-based determinations of the 
r~| ■ leading order hadronic contribution to a^. 

D 



PACS numbers: 13.40.Em,14.60.Ef,13.66.Bc,13.44.Dx 



^ ■ I. INTRODUCTION 



It is well known that, after the large, purely leptonic contribution, the largest of 
the remaining Standard Model contributions to the anomalous magnetic moment of the 
muon, = (g — 2) M /2, is that due to the leading order hadronic vacuum polarization, 
l a ^had- This contribution may be evaluated, in terms of experimental e + e~ — > hadrons 
cross-section data, using the dispersion integral representation |l[ 

lo «Im(0) r ,K(s 
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where the form of K(s) is well-known |l( and R[s) is the ratio of the "bare" e + e _ — > 
hadrons cross-section to that for e + e~ — > p + p~ |2|. Since the isovector part of the elec- 
tromagnetic (EM) spectral function is related by CVC to the charged current isovector 
vector spectral function, which can be obtained from the invariant mass distribution of 
states with zero net strangeness in the decay in r~ — > v T + hadrons, the high preci- 
sion hadronic r decay data of Refs. [j, 0] can, in principle, be used to improve the 
determination of [a/jf^ pFL IsL lol- liol Eq . 

The high accuracy achieved by the current experimental determination of a M 12l 



places a significant premium on reducing the error on [a^ had i which currently dominates 
the uncertainty on the Standard Model (SM) prediction for a M (see Ref. jll| for a recent 
review) . At the desired level of precision, the r decay data can be used only after taking 
into account the small isospin-breaking (IB) corrections to the CVC relation between 
the charged and neutral current isovector spectral functions. A detailed investigation 
of possible sources of such corrections, for the numerically dominant txtx contribution, 
has been made in Refs. 0, 3|> and the resulting s-dependent IB correction factor 



incorporated into the latest r-based analyses [3, [l6| of [ffy]^- A comparison of the 
corrected, r-based spectral data with that obtained from the recent high-precision CMD- 
2 experiment 17[ , however, shows significant residual disagreement in the tttt components 
of the two versions of the isovector spectral function ^| : the two are compatible below, 
and in the vicinity of, the p peak, but differ by ~ 5 — 10% for between ~ 0.85 
and ~ 1 GeV |l5|. This discrepancy leads to incompatible determinations of [cJ/J^, the 
r-based determination lying ~ 2a higher [l5|, 3, 18| , and producing a SM prediction for 

in agreement with the experimental result, while the EM-based determination yields 
a SM prediction which differs from experiment by ~ 2.5a (see Ref. and references 
therein for more details). 

The preliminary KLOE e + e _ — > -it + it~ radiative return data supports the earlier 
EM-based determination, yielding; a value of [a^]^ compatible with that obtained us- 
ing the CMD-2 irir data [l6|. However, as has been pointed out by many authors, the 
point-by- poi nt agreement between the CMD-2 and KLOE cross-sections is less than sat- 
isfactory [l6|, the KLOE data lying higher than CMD-2 below the p peak and lower than 
CMD-2 both on the peak and above it. The structure of the weight K(s)/s is such that 
the effects of these discrepancies largely cancel in [o^]^, but the situation nonetheless 
remains unsatisfactory. 

Differences in the p° and p ± masses and widths, suggested as one possibility for re- 
solving the TiTi spectral function discrepancy appear able to reduce locally, but not 
resolve fully the discrepancy |l6| |. 

Recent developments further complicate the picture. In Ref. jlJ], QCD sum rule 
constraints on the electroproduction and r decay data were investigated. Sum rules of 
the form 

"SO 



f"° -1 f 

/ dsw(s) p(s) = / dsw(s)Il(s) 



' s th J\s\=so 

were employed, where II(s) is either the EM or charged isovector vector current correlator, 
p(s) is the corresponding spectral function, s t h is the relevant threshold, and w(s) is a 
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function analytic inside and on the contour \s\ = Sq. The OPE is employed on the 
RHS, providing the desired constraints. At the scales employed, the OPE for the vector 
current correlators is essentially entirely dominated by the dimension D = perturbative 
contribution, and hence determined by the single input, a s . This input may be taken 
from high-scale determinations of a s (Mz) which are independent of the EM and r data 
being tested. It turns out that both the normalization and so-dependence of the weighted 
spectral integrals generated from the hadronic r decay data are in excellent agreement 
with OPE expectations In contrast, the weighted EM spectral integrals (obtained 
using CMD-2 data for the 7T7t spectral component) do not agree with OPE expectations, 
having (i) normalizations which are ~ 2a low, and (ii) slopes with respect to s which 
are ~ 2.5a low [21]. These observations suggest either a problem with the EM data, 
or the presence of non-negligible non-one-photon physics contributions to the EM cross- 
sections. In either case, the results favor determinations of [a/J^ which incorporate 
hadronic r decay data over those based on EM data alone, and a SM prediction for 
in agreement with the current BNL experimental result 0]. The recently-released SND 
e + e~ — > 7r + 7r~ cross-section results (22| are compatible with the IB-corrected r data, and 
support this conclusion. 

In light of the above unsettled situation, we revisit the question of the reliability of the 
determination of the IB corrections which must be applied to the r decay data, focussing 
on two aspects of the existing treatment. We denote the correction to [tfy]^ associated 
with these IB corrections by [Sa^]^. 

The first point concerns the uncertainty on the estimate for the contribution to [5a^]^ d 
associated with u p-u mixing" [23| (present in the EM, but not the r, spectral function). 
The most recent updates of the r-based evaluation of [a a 

] L h ad B 0, employ the IB 



corrections of Ref. |lj| (CEN). The CEN analysis is based on a version of the ChPT- 
constrained model for F n (s) developed by Guerrero and Pich |25[ (GP). The original GP 
model, which involved only the isospin conserving (IC) component of F^^s), was modified 
by CEN through the addition of an IB uj — > tttt contribution having the nominal p- 
uj mixing form. We refer to the resulting model as the GP/CEN model. Using the 
parameter values given by CEN, that part of the full IB correction associated with p-uj 
mixing becomes 

[^C m » = (3-5±0.8)xl0- 10 , (3) 

with the quoted uncertainty due essentially entirely to the 20% uncertainty on the pa- 
rameter 6 pu) , which describes the overall strength of p-uo "mixing" in the model. The 
uncertainty in Eq. (j3J) represents only a minor component of the total ±2.6 x 10~ 10 
uncertainty quoted by CEN for the full set of IB corrections ■ 



The GP/CEN mo del, however, is not the only one available for F n (s). The Gounaris- 
Sakurai (GS) model [26L t he Kuhn-Santamaria (KS) model and the hidden local 
symmetry (HLS) model |28j], for example, all predate the GP/CEN model and have been 
used extensively in the literature. The models differ in the form employed for the broad 
IC component of the e + e~ — ► 7r + 7r~ amplitude, which is given (or dominated) by the 
e + e~ — > p° — > 7r + 7r~ contribution. Implicit in the CEN error estimate is the (given the 
narrowness of the u, plausible) assumption that the value obtained for [^o, ll ] had . rnix will 
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be largely insensitive to which of the models is employed in extracting the interference 
signal. The high level of cancellation in the i^(s)/s-weighted integrals of the interfer- 
ence components of the various model cross-sections, however, makes [Sa^had-mi muc h 
more model-dependent than would be naturally anticipated. The resulting theoretical 
systematic error turns out to significantly exceed that associated with uncertainties on 
the fitted model parameters for any given model, including the GP/CEN model. This 
point is discussed in more detail in section ITD below. 

The second point concerns an IB correction not accounted for in the CEN analysis. 
In the limit that (as for the IC component) the IB component of the e + e~ — ► 7r + 7r~ 
amplitude is assumed dominated, away from threshold, by resonance contributions, three 
such contributions will, in principle, be present in the p, uj resonance region. These are 
shown in Fig. ^ where the open circles represent IC vertices and the crossed circles IB 
vertices. Jf. and Jf. are the isovector and isoscalar members of the vector current octet. 
The first two graphs represent the p-u mixing and direct IB uj — > nn decay contributions 
to the amplitude. They are small away from the uj peak region, generate contributions to 
the flavor '38' part of the EM spectral function, and combine to produce the prominent 
narrow interference shoulder in the experimental cross-section. The remaining graph 
depicts the contribution associated with the IB (isoscalar) component of the p° EM 
decay constant. Such a component of the decay constant is unavoidable in the SM. 

Because of the narrowness of the p-uj interference shoulder, the interference part of the 
cross-section can, modulo the model dependence noted above, be determined experimen- 
tally. The corresponding contribution to can thus also, with the same caveat, 
be determined experimentally. In contrast, the interference contribution associated with 
the isoscalar p° EM decay constant (which also belongs to the flavor '38' part of the 
cross-section) is identical in shape to the dominant, broad, IC flavor '33' contribution 
and not (even in principle) extractable experimentally. The corresponding contribution 
to [Sa^],,, which is certainly present at some level in the SM, has not, to our knowledge, 
been investigated in the literature, and certainly is not included in the treatment of IB 
corrections employed by CEN, for the reason explained below. It is, in fact, analogous to 
the p-uj "mixing" contribution, which was also not present in the GP model approach j 1 31 ] - 
and hence had to be added by hand to the GP model expression by CEN |l4j ]. 

The reason the three IB resonance contributions to F n (s) shown in Fig. are not 
incorporated in the GP model framework is as follows. The GP model is constructed 
by implementing the constraints of unitarity, analyticity and short-distance QCD, and 
requiring that the model expression for F n (s) match properly onto the known next-to- 
leading order (NLO) ChPT expression at low energy. This last constraint is realized using 
the resonance chiral effective theory approach, in which low-energy resonance effects 
appear through contributions proportional to the NLO low energy constants (LEC's), 
L r k , of Gasser and Leutwyler 3(|. It is, however, straightforward to demonstrate that, 
at NLO, (7r + 7r~| J^|0) receives contributions only from loops, and not from the NLO 
LEC's |3l|. As a result, none of the resonance-induced IB effects depicted in the figure 
are incorporated in the GP expression for F w (s). Although, numerically, the IB loop 
effects are tiny near threshold [3jj, the obvious experimental interference shoulder in the 
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FIG. 1: Isospin-breaking resonance contributions to e + e — * tt + tt 

p-uj region shows that this does not remain the case at higher energies. In order to include 
p-uj interference, it was thus necessary for CEN to add a p-u "mixing" contribution to 
the GP model by hand. The broad IB p contribution is similarly absent in the GP model 
approach, and would also have to be added by hand. In Section 11111 we investigate sum 
rule constraints on this contribution, and, in addition, use the size of the analogous effect 
in the pseudoscalar sector, as evaluated at NLO in the chiral expansion, to obtain some 
guidance as to what the natural scale of the effect might be. 
A brief summary, and our conclusions, is given in Section IIVI 



II. MODEL DEPENDENCE OF \5a u ][°, 

L A 1 J nad:inix 



The pion form factor, in the GS model, is given by |2( 

^ GS) (*) = (BWf s \s) [l + S^PM)] + (3BWf\s)^ (4) 



6 



where 



BW 



PJs) 



(GS), 



s - im„,r„ 



m\ 



(mi 



s + f(s, my, Ty) — im v Ty(s, my, Ty)) 



(5) 



with 



d(m v ) 
f(s,m v ,T v ) 
His) 



3m 2 



{7ipl(m 2 v )) 



hi 



(m v + 2p w (my)) 
2m n 



+ 



m v 



m%m v 



(2upjm 2 y)) (npl{m 2 v )) 



pi(s)[H(s) - H(m v )] + (m v - s)p$(m 



v) 



dH 

ds 



[m, 



£n ( Vj+ %pM 
2m^ 



TTWS 



(6) 



where p n (s) = </| — m 2 is the pion CM momentum for squared invariant mass s, 
Ty(s, my, Ty) is the standard s-dependent width for vector meson V implied by p-wave 
phase space, and Ty = Ty{m v , my, Ty). 
Similarly, for the KS model, one has j2?| 



F (KS ){i 



1 + /5 + 7 



(7) 



with 



Pv(s) 



rn 



v 
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v 



s — imyT v (s, m v , T\ 



The HLS model 28], as implemented by CMD-2, has the form 

5PM) 



Fi HLS \s) 



1 _ aHLS cihls ( Pp[s) ( I 



1 + 5 



(9) 



with anLS a constant. The model provides a good quality fit to the data below 1 GeV 
despite having no explicit p' contribution. It also turns out to rep roduce the correct final 
state 7T7T phases after the model parameters have been fitted |29| . 

For all of the GS, KS and HLS models, the constant 5, which parametrizes the strength 
of the narrow IB amplitude, is taken to be complex. A non-zero phase is, in general, 
unavoidable in the presence of an IB direct u — >• nn decay contribution |23| . 

The GP model for the pion form factor is given by [2J 



Fi GP \s) = PJs) exp 



967T 2 f 2 



ReL 



m; m; 



mj 



+ - ReL 



m 
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m 



K 
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m 2 p 



(10) 
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where 

' m 2 m 2 \ n fm 2 \ 8m 2 5 . s . f/3(s) + ll 

— ,— )=en — + +/3 S 3 ^ g 11 

s 77i*/ \ m p J s 3 Lp(s) — 1 

with /3(s) = a/1 — Aw? I s and the s-dependent width, T p (s,m p ,T p ) appearing in P p (s) 
replaced by the resonance chiral effective theory expression 

r p (a) = ^ (0(s - Ami) PA*) 3 + \o{s - 4rr£) /M*) 3 ) . (12) 

Some IB effects are incorporated into F^ GP \s) if one evaluates the phase space factors 
in the s-dependent width using the physical charged 7r and K masses. In Ref. ^H, a 
small rescaling of the coefficient appearing on the RHS of Eq. ()12j) is allowed in order to 
account for the ~ 1.5 MeV contribution of im'j decays to the total width of the p |lj, 32]. 
The CEN modification of F^ GP \s), designed to incorporate the p-uj mixing contribution 
not included in the original GP model, then has the form 

Fr /cEN \s) = rr\s) - p P (s) (^) (jl^ p^ S ) . (13) 

The parameter 9 puj was assumed real by CEN. 

The original version of the GP/CEN model, as parametrized by CEN, predates the 
most recent, corrected version of the CMD2 data, and does not provide a good fit to it, 
producing a x 2 of 80 for the 43 CMD-2 data points. The fit quality can be improved by 
allowing m p , 8 puJ (still assumed real) and the rescaling of the resonance chiral effective 
theory width to be fit to data, but the resulting optimized fit still has a \ 2 °f 61 for 
the resulting 40 degrees of freedom. The corresponding mixing contribution, [<5<v]^. mi , c > 
is shifted only slightly from the original CEN value, from 3.5 x 10~ 10 to 3.7 x 10~ 10 . 
The fact that the p-u "mixing" signal is actually a combination of mixing and direct 
uj — > 7T7T effects, however, means that an effective representation for the combination 
of the form given by Eq. (JT3j) is not generally possible without allowing 9 puJ to have a 
non-zero phase 24]. If we extend the GP/CEN model in this way, treating the phase as a 
fourth parameter to be fit to the data, an acceptable fit, having x 2 — 41 for the remaining 
39 degrees of freedom, becomes possible. We refer to this version of the GP/CEN model 
as GP/CEN* in the table below. 

The values obtained for [&v]j. . m the various models are given in Table HI All 
results are generated using versions of the models optimized to the most recent versions of 
the CMD-2 bare cross-section data [33]. Only the GP/CEN* version of the GP/CEN 
model is included since the unmodified version does not produce an acceptable quality 
fit. 

Two things are evident from the table. First, the sensitivity of \& a v\h^ d . m i X to the model 
employed is much larger than that associated with the uncertainties on the values of the 
fitted model parameters for a given model. Second, by comparing the GP/CEN* results 
to those for the GP/CEN version of the model (quoted above), we see that allowing a 
non-zero phase for IB parameter 6 pw leads to a significant decrease in [^Qjjw.™^- (A 
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TABLE I: [<5aJ, , ■ for the various models discussed in the text fit to the most recent CMD-2 
bare e + e~ — > -kit cross-sections 



Model 


X 2 /dof 


[tfajf^ • x 10 10 

L had\mix 


GS 


36/38 


2.0 ±0.5 


HLS 


37/38 


4.0 ±0.6 


KS 


37/38 


3.9 ±0.6 


GP/CEN* 


41/39 


2.0 ±0.5 



similar effect is produced by the phase of 5 in the other models.) The origin of these 
effects is easy to understand. In the interference region, the IC amplitude is ~ B p (s), 
where B p (s) is the p(770) Breit-Wigner-like form in the given model. Writing the IB 
amplitude in the u p-u mixing" form, generically B p (s)5P UJ , with 5 = \5\e l( ^, the flavor '38' 
component of the EM cross-section is then given approximately by the expression 



\B P (s)\ 2 



' " (cos(0) (ml - s) - muTusin(4>)) 



[{ml - sf + mlTl] 



(14) 



Since the coefficient multiplying cos(<f)) in the square bracket of Eq. ()14j) is antisymmetric 
about s = ml, the corresponding contribution to [$ a n\ had . rnix vanishes in the limit that 
one neglects the variation of \B p (s)\ 2 K(s)/s over the to region. Since both \B p (s)\ 2 
and K(s)/s are decreasing functions in this region, a small residual positive contribution 
remains. The coefficient of sin(<f)), in contrast, is symmetric, so no analogous cancellation 
is present in the corresponding contribution to [^ a fi)had-mixi he -i ^ ne s ^n(0) integral is 
strongly enhanced relative to the cos(0) integral. It is the strong cancellation in the 
cos(<j)) integral, combined with small differences in the s-dependence of B p (s) in the 
different models, which accounts for the significant model dependence in the results for 
[^ a fj\had-mix- Since fits to the data favor small positive <\> for all the models considered 
here, the sin(4>) contribution to [Sa^ d . ■ is negative. The relative enhancement of 
the sin(4>) integral means that the cancellation against the cos(0) contribution, which is 
absent if one sets to zero from the outset, can be quite significant, even for relatively 
small 0. 

Table HU shows the impact of the choice of input data set on [<5Gt M ]^. mix , giving the 
central values corresponding to optimized fits of each model to the bare CMD-2 and SND 
cross-sections. In the case of the KLOE data, the optimized fits for all four models have 
X 2 /dof > 2 35]; the corresponding [Sa^had-mix resu ^s have therefore been omitted from 
the table. A x 2 /dof > 2 is also obtained for the optimized fit of the GP/CEN* model to 
the SND data; the corresponding entry in the table, though included for completeness, 
has been enclosed in parentheses to remind the reader of this fact 0] . From the table we 
see that the variation in {fi a p]had m ix V8 ^ nes among the different models, for a fixed input 
data set, is significantly larger than the variation of the results for a given model over 
the different input data sets. It is thus the theoretical systematic error associated with 
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TABLE II: Central values of [&Z;,,1i. j ; , in units of 10 10 , for the optimized fits of the various 
models to the bare CMD-2 [ljj and SND |22| e + e~ — > 7T7r cross-sections. 



Model 


CMD-2 


SND 


GS 


2.0 


2.2 


HLS 


4.0 


4.5 


KS 


3.9 


4.3 


GP/CEN* 


2.0 


(1.6) 



TABLE III: Central values for [Sap] had . mix for earlier model fits in the literature 



Model Reference 


x 2 /dof 


MS***, x 10 10 


GS ALEPH97 


59/38 


4.3 


KS ALEPH97 


75/38 


6.3 


GS DAVIER03 


65/38 


2.5 


GP/CEN CEN02 


80/42 


3.5 



choice of model used in separating the IC and IB components of the amplitude which 
dominates the uncertainty in the determination of [<5a M ] jLw- ; ■ 

For comparison, and to further illustrate the sensitivity of the mixing contribution 
to small changes in the data and the resulting model fits, results corresponding to some 
older fits from the literature are given in Table ITTT1 'CEN02' labels the original GP/CEN 
result jlH (with no phase for the parameter 9^), 'ALEPH97' the results corresponding 
to the "combined" (r plus electroproduction) GS and KS fits of Ref. jij and 'DAVIER03' 
the results corresponding to the similarly "combined" GS fit (Table 4) of Ref. [43j|. Details 
of the fit procedures, and data sets employed may be found in the original references. 
The total x 2 of the fits relative to the 2003 CMD-2 data are also given. The quoted 
results correspond in each case to the central values of the fit parameters for the models. 
Typically, differences between the old and new values of [<fo M ] /wrf . mia . f° r a given model 
are much larger than one might have anticipated, given the relatively small changes in 
the both data and fitted parameter values. This sensitivity is again a reflection of the 
strong cancellation in the integral of the product of IC and IB amplitudes. 

We conclude that the inability to separate the IC and IB components of the e + e~ — > tttt 
amplitudes in a model independent manner leads to a significant uncertainty in the 
evaluation of [Sa^]^ d . mix . How one assesses this uncertainty depends on one's attitude 
to the various models. One stance might be to argue that the GS and GP/CEN models, 
which explicitly incorporate the constraints of unitarity and analyticity, are to be favored 
in deciding on a central value. The variation of the results across the different models 
would then serve as a measure of the residual uncertainty. Alternately, since all of the 
models are, to greater or lesser extent, phenomeno logical, and, at least for the CMD-2 
data, yield comparable quality fits, one could instead average the results to arrive at 
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a central value, and assign an error large enough to incorporate the highest and lowest 
values allowed by the errors associated with those on the fitted parameters for the various 
models. The first stance would yield 

WEinU. = (2.0±2.9)xl0- 10 , (15) 

the second 

= (3.1±1.8)xKT 10 . (16) 
A smaller central value, and significantly larger uncertainty, would result if one ignored 



the poor quality of the model fits and also took the KLOE-based results into account (37 



III. THE BROAD IB p CONTRIBUTION 

Although the broad IB p contribution to the experimental e + e~ — > 7t + tt~ cross-section 
has the same shape as the dominant IC contribution, and hence cannot be separated from 
it experimentally, it can, in principle, be determined theoretically through a QCD sum 
rule analysis of the IB vector current correlator, II 38 (q 2 ), defined by 

• J d'xe^ (0|T (4(x)4(0)) |0> = (q,q v - q 2 g,„) II 3 V) . (17) 

A non-zero (IB) coupling of the p to J 8 will produce a broad p contribution to the 
spectral function of II 38 , p 38 (s), whose strength, X p , is proportional to the product 
of the IC isovector and IB isoscalar p° decay constants. This product, together with 
analogous IB products, Xy, for the other vector meson resonances, can, in principle, be 
determined by matching the appropriate weighted integrals of p 38 to the corresponding 
OPE expressions. 

In Ref. |38| such an analysis was performed using two different families of "pinch- 
weighted" finite energy sum rules (pFESR's) j^. The details of the analysis may be 
found in Ref. (38[, and will not be repeated here. The following point is, however, worth 
noting. Because the numerically dominant term on the OPE side of the various sum 
rules is that with dimension D — 4, the contribution from the VEV's of D = 6 four- 
quark operators, cqOq, is not expected to be negligible. The fact that the VEV's for such 
operators are typically not well known empirically would normally present a problem for 
the sum rule analysis. It turns out that the dependence of the Xy on cqOq is different 
for the two different pFESR families, allowing, not only the Xy, but also cqOq, to be 
determined from the combined analysis. The values of CqOq which make the different Xy 
consistent turn out to agree at the ~ 1% level, providing strong support for the reliability 
of the analysis. 

Unfortunately, in the analysis of Ref. [3^, the p-u "mixing" contribution to p 38 (s) 
implied by the observed interference shoulder in the EM cross-section was not input 
separately on the spectral side of the sum rules employed. As a result, the output X p 
contains contributions from all three of the IB processes shown in Figure In order 
to separate out the experimentally inaccessible contribution associated with the third 
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of those processes, we have redone the analysis of Ref. |38[, this time inputting the 
"interference" component of p 38 , as determined from the CMD-2 experimental data in 
the interference region. This input, as noted above, depends to some extent on the choice 
of model for the p contribution to the IC component of the amplitude. We have then used 
the various pFESR's to solve for the residual broad p contribution, which yields directly 
the contribution to the flavor '38' part of the EM spectral function associated with the 
third graph in Fig.0 With current experimental errors, an accurate determination of this 
"direct" contribution, and hence of the associated contribution to [&x M ]f,, [Sa^had-direcv 
turns out to be impossible. We find central values of a natural scale (see below), but with 
errors, induced by the uncertainy in the integrated (model-dependent) p-u "interference" 
term, much larger than these central values. Even more unfortunately, versions of the 
analysis using errors scaled down artificially by hand suggest that the reduction needed 
to allow even just a reliable determination of the sign of the direct effect are unlikely to 
be reachable in the foreseeable future. 

In view of the weakness of the constraints arising from the sum rule analysis, we 
turn to the pseudoscalar sector, and study the size of analogous effects in the IB decay 
constant, ft, of the 7r°, defined by 

(0|4|7T°> = if,e iq , = ifX- (18) 

At leading order in the chiral expansion, the IB parameter e± = ft/ f-K is equal to 9 = 
V^^d — m u ) / A(m s — m), where m = (m^+m u )/2, and is due entirely to IB mixing on the 
external leg, induced by the non-zero value of md — m u . At NLO, e\ receives contributions 
both from mixing and from IB in the low energy representation of the axial current A^. 
The full NLO expression for e± is given in Ref. [3(j, while the ingredients necessary for 
separating the mixing and vertex contributions may be found in Ref. At NLO the 
separate mixing and vertex contributions ctrc , cts expected on physical grounds, scale 
dependent. The (scale-dependent) "direct" (vertex) contribution is given by 



ft 



direct 



16m 2 K + 4ml\ + g m| m\ £n(m 2 K / p 1 



ft ) ^' ft ^' ^ 2 ft 



+ l ! il^i £ J [1+M<//i2)] ' (19) 

where p is the ChPT renormalization scale, and L r 45 (p) are the renormalized NLO LEC's, 
evaluated at scale p [30j | . 

Taking, for illustration, the central values for LJ S from the main fit of Ref. [4lj], one 
finds that the leading (NLO) contribution to ft varies from 0.58% to 0.46% of f w as a 
varies from to 1 GeV. Using the slightly modified fit denoted "fit D" in Ref. |42J, 
which produces good values for the threshold parameters for tttt and ttK scattering, the 
ratio runs from 0.60% to 0.71% over the same range of p. We thus conclude that the 
NLO "direct" contribution to the IB decay constant ratio is ~ 1/2% for typical hadronic 
scale choices. A similar value for the corresponding IB ratio of p° decay constants would 
produce a contribution \& a v\^ d . diTect — (2 — 3) x 10~ 10 . The central value obtained from 
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the sum rule analysis turns out to be comparable to, or less than, this natural size for all 
four models considered here. The upper bound implied by the errors, however, is much 
larger, preventing us from using the sum rule constraints in a meaningful way. 

Note that the presence of a small IB component in the p° EM decay constant would 
have an impact on the values of m p and T p obtained by fitting the various phenomeno- 
logical models to data. At the natural scale (given above) for this effect, however, the 
(model-dependent) shift in m p would be ~ 0.2 MeV or less, and that in T p 0.4 MeV or 
less. The effect, while contributing to the apparent difference in charged and neutral 
p masses and widths, can thus account for at best a modest fraction of the differences 
obtained after fitting the models independently to the r decay and electroproduction 
data. 



IV. CONCLUSIONS 

We have shown that the model dependence encountered in separating the IB from 
the IC component of the e + e~ — ► 7r + 7r~ cross-section leads to a theoretical systematic 
uncertainty on the p-to "mixing" contribution to [Sa p ]^ d which is several times that 
associated with the fit parameter uncertainties for any given model. There is also a 
potentially non-negligible IB contribution, associated with the direct IB coupling of the 
isoscalar part of the EM current to the p°, which cannot, even in principle, be determined 
experimentally. At present, we are able only to make a rough "natural size" estimate 
for the magnitude of this contribution. Both of these effects would need to be taken 
into account when estimating the uncertainty in the IB correction required in order to 
incorporate the r~ — > v T ix~ix® spectral data into the evaluation of [o/J/Ld- The results 
make it clear that evaluating the "mixing" component of the IB correction using only a 
single model in the analysis of the cross-section will lead to a significant underestimate 
in the uncertainties. 

In view of the results of the tests involving independent high-scale OPE constraints 
reported in Ref. [2l[, and the compatibility of the SND and r nn data sets, we base our 
final results on the SND data set. (The CMD-2 and KLOE values, quoted above, allow 
alternate choices to be made.) Choosing the GS model as the favored case would lead 
to a decrease of 1.3 x 10 -10 in [<5ajf°, . , relative to the standard CEN value. Such a 

shift would lower the r-based prediction for [5a M ], , in the SM, slightly increasing the 
difference between the SM prediction and the central experimental value, but leaving 
them compatible at the la level. For a given model, the results for [m w . mja . obtained 
using either the CMD-2 or SND data are actually in good agreement. The decrease in 
central value relative to CEN would thus also reduce the discrepancy between the CMD-2 
and r-based determinations of the the tttt contribution to [<5<^]^ by ~ 10%. The HLS 
model, which produces an optimized fit to the SND data of comparable quality to that 
of the GS model, in contrast, shifts the mixing contribution up by 1.0 x 10~ 10 , relative 
to the CEN value. Since the GS model has the constraints of analyticity and unitarity 
explictly built into it, while the HLS model does not, we have favored the GS model in 
our discussion of the central value above. It is, however, important to bear in mind that, 
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in terms of quality of fit to the SND data, the models cannot be distinguished, and hence 
that a significant uncertainty, associated with the model dependence, must be attached 
to any particular chosen central value. The range of values allowed by the CMD-2 data is 
fully contained within that allowed by the SND data, but the latter range would need to 
be significantly extended, including to negative values of [<5<v]^. mix , if the lower-quality 
fits to the KLOE data were also taken into account. 

Finally, taking into account the uncertainties already identified by CEN, and adding to 
these both the increased uncertainty on [^a^jfj. . and a possible direct IB isoscalar p° 
EM coupling contribution of the "natural size" discussed above, we arrive at a combined 
uncertainty for the IB correction one must apply in order to use the t~ — > z/ T 7r~7r° data 
in computing [(foyjf . which is ~ 4 x 1CT 10 . It appears unlikely that this uncertainty can 
be significantly reduced. Should the new BNL experimental proposal 0] be approved, 
the uncertainty on the IB corrections would thus exceed those on the experimental de- 
termination of a M , seriously limiting our ability to make use of hadronic r decay data in 
determining the SM prediction for [5a M ]^. 
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